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Abstract — Linear minimum mean square error (MMSE) detec- 
tor has been shown to alleviate the noise amplification problem, 
resulting in the conventional zero-forcing (ZF) detector. In this 
paper, we analyze the performance improvement by the MMSE 
detector in terms of the condition number of its filtering matrix, 
and in terms of the post-precessing signal to noise ratio (SNR) 
improvement. To this end, we derive explicit formulas for 
the condition numbers of the filtering matrices and the post- 
processing SNRs. Analytical and simulation results demonstrate 
that the improvement achieved by the MMSE detector over the 
ZF detector is not only dependent on the noise variance and the 
condition number of the channel matrix, but also on how close 
the smallest singular values are to the noise variance. 

I. Introduction 

Multiple-input multiple-output spatial multiplexing (MLMO- 
SM) techniques are considered as the key technology to in- 
crease the channel capacity in mobile communication systems 
beyond third-generation (3G) (TJ. In MIMO-SM, signals are 
transmitted simultaneously from different spatial elements, 
i.e., antennas, using the same time-frequency resources Q. 
Therefore, the channel capacity is increased without requiring 
any additional spectral resources, which are not only expensive 
but also scarce. Indeed, the attained capacity by the MIMO 
multiplexing techniques depends crucially on the detection al- 
gorithm which is employed at the receiver side to demultiplex 
the transmitted signals. 

Maximum-likelihood (ML), which employs brute-force 
search, is the optimum detection technique for the MIMO mul- 
tiplexing systems |3|. Nonetheless, its exponential complexity 
is inapplicable in computational complexity and latency lim- 
ited mobile communication systems. 

Several suboptimal detection schemes were proposed in the 
literature [4|-[7|. It is clearly shown in several works that 
taking the noise into consideration increases the detection 
efficiency (see for instance [8| and references therein). This 
improvement in the efficiency can be translated into reduc- 
tion in the computational complexity, error enhancement, or 
both of them. Nonetheless, this improvement achieved by the 
minimum-mean square error (MMSE) criterion is referred to 
the value of the noise variance, or more accurately to the 
signal to noise ratio (SNR). That is; when the noise variance is 
small, the performances of the ZF and MMSE linear detectors 
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coincide J5]. We show that this reasoning, which is frequently 
used in the literature, is inaccurate. 

Contributions. In this paper, we investigate the improve- 
ment achieved by the linear detector using the MMSE criterion 
over that using the ZF criterion. The performance improvement 
is given as the ratio between the condition numbers of the 
filtering matrices. We show that the achieved improvement 
by the MMSE criterion does not only depend on the noise 
variance, but also on how close the noise variance is to the 
singular values of the channel matrix. Also, we derive the post- 
processing SNRs of the linear ZF and MMSE detectors and 
show the superiority of the detection algorithm when noise is 
considered in constructing the filtering matrix. 

The rest of this paper is as follows. In Section II, we 
introduce the system model and give a brief review of linear 
detection schemes. In Section III, we analyze the condition 
numbers of the filtering matrices of the ZF and MMSE detec- 
tors. Analytical and simulation results are shown in Section 
IV, and conclusions are drawn in Section V. 

We define now some variables and notations used through- 
out the paper. The notation Tr(-) denotes the trace of a matrix, 
E(-) denotes the expectation, and (-) H is the Hermitian trans- 
pose. A;(H) and <7;(H) are the i-th eigenvalue and singular 
value of the matrix H, respectively. 

II. System Model and Review of the Linear 
Detection Techniques 

We consider a MIMO spatial multiplexing system with nr 
transmit and tir receive antennas. For the sake of simplicity, 
we assume that N = tit = ur. The received vector r € 
is given by: 

r = Hx + n, (1) 

where n is the Gaussian noise vector with E 
and H is the channel matrix whose element hij is the zero- 
mean unit-variance complex Gaussian random variable. Also, 
x is the transmitted data vector whose elements are drawn 
independently from a quadratic amplitude modulation (QAM) 
set. 

Linear detection algorithms are the simplest approaches 
used to solve (Q}. The main idea behind the linear detectors 
is to treat the received vector by a filtering matrix W, which 
is constructed using a performance-based criterion. The zero- 
forcing and minimum-mean square error criteria are used in 
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Fig. 1. MIMO spatial multiplexing with linear receiver. 



the linear ZF and MMSE detection algorithms, respectively. 
Hence, the filtering matrices for the ZF and MMSE detectors 
are given as follows: 
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where a\ is the noise variance and Ijy is the N x N 
identity matrix. Fig. Q] depicts the MIMO multiplexing system 
with linear receiver. In the sequel, the channel matrix H is 
considered to be perfectly known at the receiver. 

III. Achieved Improvement by the MMSE Detector 

A. Improvement in the Conditionality of the Filtering Matrix 

In this section, we evaluate the improvement achieved by 
the linear MMSE over the linear ZF algorithm in terms of 
the condition number of their corresponding filtering matrices. 
Note that the condition number is used to demonstrate the 
accuracy of the recovered signal, and to show the system's 
immunity against channel estimation errors. The condition 
number of an M x TV matrix A is defined as the ratio between 
its maximal and minimal singular values; that is, 



cond(A) = ||A|| • II A" 1 1| = 



t(A) 



i(A)' 



(4) 



where ||A|| is the Euclidean norm of the matrix A, which 
equals the maximum singular value of A. To proceed in the 
analysis, we introduce the following Inequality and Lemma. 

Inequality III.l. Weyl's Inequalities (Eigenvalues) iflOl : Let 

£ and A be N x N Hermitian matrices, Ai(£) > A2(£) > 
• • • > Aat(S) > 0, Ai(A) > A 2 (A) > • • > Ajv(A) > 0, and 
Ai(r) > A2(r) > • • • > AAr(r) > are the eigenvalues for 
£, A, and T = (£ + A), respectively. Then, 
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Lemma III.2. Let A be an N x N Gram matrix, i.e., A = 
B ff B for any nonsingular matrix B. Then, \(A) > and 
Ai(A) is also a singular value. 

Proof: By definition, ct ( (H) = A, 1/2 (H H H), which im- 
plies that cr?(H) = A^H^H). Notice that the singular 
value decomposition (SVD) of (H H H) = UE^V^VEU^ = 



U(S if E)U H , which is equivalent to the eigenvalue decom- 
position with orthonormal eigenvectors. Therefore, erf (H) = 
(Ji(H ff H) = Ai(H H H). ■ 
From Weyl's inequalities for the eigenvalues and Lemma 
dlinil, it follows that: 
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Now, we define the Gram matrices £ = K H K and A = A H A, 
where A = <7 n lAr. Then, dividing both sides of (0 by the 
corresponding sides of © results in: 
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Lemma III.3. For any nonsingular N x N matrix A, 
cond(A) = cond(A _1 ). 

Proof: The proof of this lemma follows from the fact that 
o~i (H — 1 ) = l/crjv(H) and ^(H^ 1 ) = l/cri(H). ■ 
Then, the ratio between the condition number of the filtering 
matrices of the linear MMSE and ZF detectors is approximated 
as follows: 
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(10) 



l + al/a^ft) 

Based on ( flOb , we make the following remarks: 

1) When the channel matrix H is orthogonal, i.e., u\ (H) = 
<7jv(H), the condition number of W z f equals that of 

Wmjusc • 

2) The improvement by the MMSE detector does not 
depend only on the noise variance cr^, but also on how 
close the noise variance is to the minimum singular value 
<7jy-(H). Also, it depends on the conditionality of the 
channel matrix H. 



B. Improvement in the Post-processing SNR 

In the case of ZF detector, the post-processing signal vector 
is given as follows: 



H r 



H n 



(11) 



The post-processing SNR of the ZF detector is then given by: 
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Let (HH H ) = V£V H be the singular value decomposition of 
the Gramian matrix (HH H ), then (fT2l can be simplified as 
follows: 



SNR zf = 



N 



where the expectation in the denominator is only taken with 
respect to the noise vector n. 

Based on the analysis carried out in ifTTl for the case 
of multi-user precoding, we conducted the analysis for the 
linear MMSE detector. The post-processing SNR is given as 
following: 
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Notice that for v\ = 0, SNR mmsc = SNR zf . 

IV. Simulation Results and Discussions 

In this section, we consider the elements of the channel 
matrix H to be independent and identically distributed (i.i.d.) 
and to follow a zero-mean unit-variance complex Gaussian dis- 
tribution. Each realization of the channel matrix is normalized 
such that 



where hi 



N N N 
i— 1 j — 1 i—1 

is the row-column element of the channel matrix H. 



Therefore, the SNR at each receive antenna equals Also, 
quadrature phase shift keying (QPSK) modulation is used. 

Fig. |2] depicts the gain attained by the MMSE detector over 
the ZF detector in terms of the post-processing SNR. As N 
increases, the gain achieved by the MMSE detector increases, 
attaining about 15dB at a receive SNR of OdB. Indeed, it is 
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Fig. 2. Achieved gain by the linear MMSE over the ZF detector in terms of 
the post-processing signal to noise ratio. 

TABLE I 

Minimum singular value and condition number of the channel 
matrix for several configurations. values are averaged over 
1 0,000 independent channel realizations . 



N 
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E(o"jv(H)) 
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0.257 


0.224 


0.198 


E(cond(H)) 


4.27 


10.82 


24.32 


37.82 


49.85 


64.95 



the post-processing SNR of the ZF detector that deteriorates 
due to the vanishing of the minimum singular value of the 
channel matrix as N increases. 

Table U depicts the expectations of the minimum singular 
value and the condition number of the channel matrix for 
several N. The given values are the average over 10,000 
independent realizations of H. As N increases, we remark 
that minimum singular value decreases, and consequently 
the condition number increases. This clearly explains the 
poor performance of the linear ZF detector, where the post- 
processing noise variance depends directly on the inverse of 
the singular values. These results coincide with the following 
theorem (see Thm. 2.36 in lfl3l and Thm. 5.1 in ifPfln . 

Theorem IV. 1. Consider the N x N standard complex Gaus- 
sian matrix H, i.e., hij ~ CN(0, 1). The minimum singular 
value of H satisfies 



N 



lim P[Na N (H) > x] = c" 



-x 2 /2 



(19) 



where P[x] is the probability of x. 

Fig. [3] shows the cumulative distribution function (cdf) of 
the minimum singular value crjv(H) for several N values. As 
expected from Table U for an N x N matrix A and an M x M 
matrix B with M > N, the following is satisfied 

P[(T M (B) <x}> P[a N (A) < x}. (20) 
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Fig. 3. Cumulative distribution function (cdf) of the minimum singular value 
of the channel matrix in several system configurations. 
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Fig. 5. Effect of the minimum singular value of the channel matrix on the 
improvement by the MMSE detector in the low SNR regime. 
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Fig. 4. Effect of the minimum singular value of the channel matrix on the 
performance of the linear detectors in 4 X 4 MIMO multiplexing system. 
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Fig. 6. Effect of the minimum singular value of the channel matrix on the 
improvement by the MMSE detector in the high SNR regime. 



This gives the reason behind the degradation in the perfor- 
mance of the ZF detector as N increases. 

Fig. |4] depicts the bit error rate (BER) performance of the 
linear MMSE and ZF detectors in 4 x 4 MIMO multiplexing 
system. Notice that the power constraint in ( TT~8T > is still valid for 
all the obtained simulations. When no constraint is imposed 
on the minimum singular value, i.e., a m i n > , the MMSE 
outperforms the ZF detector by about 4.2dB for all the 
simulated SNR range. Notice that even at SNR of 50dB, viz., 
0-2 — 2 x 1CP 5 , the improvement by the MMSE detector is 
kept fixed. This indicates that the improvement by the MMSE 



in the error rate does not only depend on the plain value of 
the noise variance. 

To show the relationship between the minimum singular 
value and the improvement achieved by the MMSE detector, 
we assign a lower bound on the minimum singular value 
of H. We remarked that the BER performances of the ZF 
and MMSE detectors coincide in the case of (<7 m i n = 0.3), 
((T min = 0.2), and (<7 min = 0.1) at SNR of 25dB, 28dB, and 
33dB, respectively. This is equivalent to cr^Ja^ equals 14, 
12.5, and 10, respectively. We don't have a good explanation of 



these results, and leave it for further studies. This clearly shows 
that the improvement by the linear MMSE detector depends 
on the ratio between the minimum singular value of H and 
the noise variance, rather than only on the plain value of the 
noise variance. 

Fig. [5] shows the ratio between the condition numbers of 
the filtering matrices of the MMSE and ZF detectors for 
1/ofj =10dB and cond(H) = 15. The "Approximation" is 
obtained via our derived formula in (TToT >. We remark that 
our approximation is accurate except for small values of 
(Tmin, where the error is tolerable. We also notice that the 
conditionality of H is highly improved at low <J m i n , where the 
filtering matrix tends to be orthogonal (cond(W mmsc ) rj 1.5) 
at cr m i n = 0.1. The improvement by the MMSE regularization 
becomes negligible as <7 m i n increases. 

Fig. [6] shows the ratio between the condition numbers of 
the filtering matrices of the MMSE and ZF detectors for 
l/al =30dB and cond(H) = 15. At this high SNR value, 
the improvement by the regularization imposed by the MMSE 
detector is also possible when er m j n is small. This explains the 
improvement in the BER for (cr m ; n > 0), as shown in Fig. |U 

V. Conclusions 

In this paper, we analyzed the performance of the linear ZF 
and MMSE detectors. Based on Weyl's Inequalities, we gave 
an accurate approximation of the ratio between the filtering 
matrices of the MMSE and ZF detectors. We revealed that the 
improvement in the conditionality of the channel matrix by 
the MMSE detector is dependent on the relationship between 
the noise variance and minimum singular values. Also, explicit 
formulas for the post-processing SNR for the ZF and MMSE 
detectors are derived, where we showed that the gain in the 
SNR by the MMSE detector increases with the increase in 
the number of transmit antennas. This is analytically argued 
by referring to the vanishing nature of the minimum singular 
value of the channel matrix as N increases. Finally, simula- 
tions were carried out, where obtained results coincided with 
our analytical conjectures and approximate formulas. 
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